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ABSTRACT 

The theory of reheating of the Universe after inflation is developed. We have found 
that typically at the first stage of reheating the classical inflaton field <fi rapidly decays 
into 0-particles or into other bosons due to a broad parametric resonance. Then these 
bosons decay into other particles, which eventually become thermalized. Complete 
reheating is possible only in those theories where a single particle <fi can decay into 
other particles. This imposes strong constraints on the structure of inflationary models, 
and implies that the inflaton field can be a dark matter candidate. 
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1. The theory of reheating of the Universe after inflation is the most important application 
of the quantum theory of particle creation, since almost all matter constituting the Universe at 
the subsequent radiation-dominated stage was created during this process [T|]. At the stage of 
inflation all energy was concentrated in a classical slowly moving inflaton field <fi. Soon after the 
end of inflation this field began to oscillate near the minimum of its effective potential. Gradually 
it produced many elementary particles, they interacted with each other and came to a state of 
thermal equilibrium with some temperature T r , which was called the reheating temperature. 

An elementary theory of reheating was first developed in for the new inflationary scenario. 
Independently a theory of reheating in the R 2 inflation was constructed in ||. Various aspects of this 
theory were further elaborated by many authors, see e.g. M. Still, a general scenario of reheating 
was absent. In particular, reheating in the chaotic inflation theory remained almost unexplored. 
The present paper is a short account of our investigation of this question 0. We have found that 
the process of reheating typically consists of three different stages. At the first stage, which cannot 
be described by the elementary theory of reheating, the classical coherently oscillating inflaton 
field (f) decays into massive bosons (in particular, into 0-particles) due to parametric resonance. In 
many models the resonance is very broad, and the process occurs extremely rapidly (explosively). 
Because of the Pauli exclusion principle, there is no explosive creation of fermions. To distinguish 
this stage from the stage of particle decay and thermalization, we will call it pre-heating. Bosons 
produced at that stage are far away from thermal equilibrium and typically have enormously large 
occupation numbers. The second stage is the decay of previously produced particles. This stage 
typically can be described by methods developed in 0. However, these methods should be applied 
not to the decay of the original homogeneous inflaton field, but to the decay of particles and fields 
produced at the stage of explosive reheating. This considerably changes many features of the 
process, including the final value of the reheating temperature. The third stage is the stage of 
thermalization, which can be described by standard methods, see e.g. fl], @; we will not consider 
it here. Sometimes this stage may occur simultaneously with the second one. In our investigation 
we have used the formalism of the time-dependent Bogoliubov transformations to find the density 
of created particles, n^{t). A detailed description of this theory will be given in here we will 
outline our main conclusions using a simple semiclassical approach. 

2. We will consider a simple chaotic inflation scenario describing the classical inflaton scalar 
field <p with the effective potential V((f>) = ±.\m 2 b cf) 2 + |0 4 . Minus sign corresponds to spontaneous 
symmetry breaking — > <fi + a with generation of a classical scalar field o = The field <fi 

after inflation may decay into bosons x an d fermions ip due to the interaction terms — \g 2 (t> 2 X 2 an d 
—h'tpipcf). Here A, g and h are small coupling constants. In case of spontaneous symmetry breaking, 
the term — ^g 2 <p 2 x 2 gives rise to the term —g 2 a(px 2 - We will assume for simplicity that the bare 
masses of the fields x an d ip are very small, so that one can write m x ((j)) = gcfi, m^((f>) = \h<f>\. 

Let us briefly recall the elementary theory of reheating At <fi > M p , we have a stage of 
inflation. This stage is supported by the friction-like term 3H<fi in the equation of motion for 
the scalar field. Here H = a/a is the Hubble parameter, a(t) is the scale factor of the Universe. 
However, with a decrease of the field <fi this term becomes less and less important, and inflation ends 
at (ft < Mp/2. After that the field begins oscillating near the minimum of V(4>) ||. The amplitude 
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of the oscillations gradually decreases because of expansion of the universe, and also because of the 
energy transfer to particles created by the oscillating field. Elementary theory of reheating is based 
on the assumption that the classical oscillating scalar field <f)(t) can be represented as a collection 
of scalar particles at rest. Then the rate of decrease of the energy of oscillations coincides with the 
decay rate of 0-particles. The rates of the processes — > XX an d 4> — » (f° r m 4> ^ 2m x , 2m^) 
are given by 

r(0-xx) = |^, r(*-v>V0 = ^- (i) 

Snm^ Sti 
Reheating completes when the rate of expansion of the universe given by the Hubble constant 
H = ~ t^ 1 becomes smaller than the total decay rate T = T((j> — > xx) + r(0 — > V'V')- The 



reheating temperature can be estimated by T r ~ 0.1 \ TM 1 



As we already mentioned, this theory can provide a qualitatively correct description of particle 
decay at the last stages of reheating. Moreover, this theory is always applicable if the inflaton field 
can decay into fermions only, with a small coupling constant h 2 <C m^/Mp. However, typically 
this theory is inapplicable to the description of the first stages of reheating, which makes the whole 
process quite different. In what follows we will develop the theory of the first stages of reheating. 
We will begin with the theory of a massive scalar field decaying into particles Xi then we consider 
the theory ^0 4 , and finally we will discuss reheating in the theories with spontaneous symmetry 
breaking. 

3. We begin with the investigation of the simplest inflationary model with the effective potential 

-^0 2 . Suppose that this field only interacts with a light scalar field x ( m x ^ m <f>) due t° the term 
— \g 2 4> 2 X 2 - The equation for quantum fluctuations of the field x with the physical momentum 
k/a(t) has the following form: 

Xk + SH Xk + + g 2 ® 2 sin 2 (m^ X k = , (2) 

where k = \f~& 2 , and $ stands for the amplitude of oscillations of the field 0. As we shall see, the 
main contribution to x - P ar ticle production is given by excitations of the field x with k/a 3> m^, 
which is much greater than H at the stage of oscillations. Therefore, in the first approximation 
we may neglect the expansion of the Universe, taking a(t) as a constant and omitting the term 
3Hxk in @ ■ Then the equation (||) describes an oscillator with a variable frequency Q 2 (t) = 
k 2 a~ 2 + g 2 § 2 sin 2 (m^t). Particle production occurs due to a nonadiabatic change of this frequency. 
Equation (|2|) can be reduced to the well-known Mathieu equation: 

xl + (A(k)-2qcos2z) X k = , (3) 

where A(k) = Jf a2 + 2q, q = z = m^t, prime denotes differentiation with respect to z. An 

<t> <t> 

important property of solutions of the equation (y) is the existence of an exponential instability 
Xk exp(fi^z) within the set of resonance bands of frequencies Ak^ labeled by an integer 
index n. This instability corresponds to exponential growth of occupation numbers of quantum 
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fluctuations n^{t) oc exp(2/^T m^t) that may be interpreted as particle production. The simplest 
way to analyse this effect is to study the stability/instability chart of the Mathieu equation, which 
is sketched in Fig. 1. White bands on this chart correspond to the regions of instability, the grey 
bands correspond to regions of stability. The curved lines inside white bands show the values of 
the instability parameter The line A = 2q shows the values of A and q for k = 0. All points 
in the white regions above this line correspond to instability for any given q. As one can see, 
near the line A = 2q there are regions in the first, the second and the higher instability bands 
where the unstable modes grow extremely rapidly, with fi^ ~ 0.2. We will show analytically in || 
that for q ^> 1 typically /i^ ~ ~ 0.175 in the instability bands along the line A = 2q, but its 

maximal value is ln ( 1+v/2 ") ~ 0.28. Creation of particles in the regime of a broad resonance (q > 1) 
with 27r/!fc = 0(1) is very different from that in the usually considered case of a narrow resonance 
(q <C 1), where 2717^ <C 1. In particular, it proceeds during a tiny part of each oscillation of the field 
<p when 1 — cos z ~ q~ l and the induced effective mass of the field x (which is determined by the 
condition m 2 = g 2 $ 2 /2) is less than m^. As a result, the number of particles grows exponentially 
within just a few oscillations of the field <fi. This leads to an extremely rapid (explosive) decay of 
the classical scalar field 4>. This regime occurs only if q > 7r _1 , i.e. for g& > m^, so that -C gM p 
is the necessary condition for it. One can show that a typical energy E of a particle produced at 
this stage is determined by equation A — 2q ~ y/q, and is given by E ~ ^ gm^M p ||. 

Creation of ^-particles leads to the two main effects: transfer of the energy from the homoge- 
neous field 4>(t) to these particles and generation of the contribution to the effective mass of the 
d> field: 

m %eff ~ m Q 2 (x 2 )ren- The last term in the latter expression quickly becomes larger 
than m 2 . One should take both these effects into account when calculating backreaction of created 
particles on the process. As a result, the stage of the broad resonance creation ends up within 
the short time t ~ mT 1 ln(m</,/^ 5 M p ), when $ 2 ~ (x 2 ) and q = j 3 -^- — becomes smaller than 

1. At this time the energy density of produced particles ~ E 2 (x 2 ) ~ gm^M p $ 2 is of the same 
order as the original energy density ~ m|M 2 of the scalar field <fi at the end of inflation. This 
gives the amplitude of oscillations at the end of the stage of the broad resonance particle creation: 
^ 2 ~ (x 2 ) ~ g~ l fn^M p M 2 . Since E 3> m^, the effective equation of state of the whole sys- 
tem becomes p m e/3. Thus, explosive creation practically eliminates a prolonged intermediate 
matter-dominated stage after the end of inflation which was thought to be characteristic to many 
inflationary models. However, this does not mean that the process of reheating has been completed. 
Instead of x-particles in the thermal equilibrium with a typical energy E ~ T ~ (mM p ) 1 '' 2 , one has 
particles with a much smaller energy ~ (gm^Mp) 1 ^ 2 , but with extremely large mean occupation 
numbers ~ g~ 2 ^> 1. 

After that the Universe expands as a{t) oc \^t, and the scalar field continues its decay in the 
regime of the narrow resonance creation q w < 1. As a result, <p decreases rather slowly, 

<p oc t _3//4 . This regime is very important because it makes the energy of the field much smaller 
than that of the x-particles. One can show that the decay finally stops when the amplitude of 
oscillations $ becomes smaller than g^ 1 m ( p ||. This happens at the moment t ~ ■m^ 1 (gM p /m ( f,) 1 ^ 3 
(in the case m < g 7 M p decay ends somewhat later, in the perturbative regime). The physical 
reason why the decay stops is rather general: decay of the particles (f) in our model occurs due to its 
interaction with another 0-particle (interaction term is quadratic in cf) and in x). When the field <fi 
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(or the number of 0-particles) becomes small, this process is inefficient. The scalar field can decay 
completely only if a single scalar 0-particle can decay into other particles, due to the processes 
~* XX or — i> ipip, see eq. ([I]). If there is no spontaneous symmetry breaking and no interactions 
with fermions in our model, such processes are impossible. 

At later stages the energy of oscillations of the inflaton field decreases as a~ 3 (t), i.e. more 
slowly than the decrease of energy of hot ultrarelativistic matter oc a~ 4 (t). Therefore, the relative 
contribution of the field 4>(t) to the total energy density of the Universe rapidly grows. This gives 
rise to an unexpected possibility that the inflaton field by itself, or other scalar fields can be cold 
dark matter candidates, even if they strongly interact with each other. However, this possibility 
requires a certain degree of fine tuning; a more immediate application of our result is that it allows 
one to rule out a wide class of inflationary models which do not contain interaction terms of the 
type of g 2 o~(j)x 2 or h(fnpip. 

4. So far we have not considered the term in the effective potential. Meanwhile this term 
leads to production of 0-particles, which in some cases appears to be the leading effect. Let us study 

HI? 1 

the 0-particle production in the theory with V{4>) = -f-4> 2 + j 4> 4 with m| <C AM 2 . In this case the 
effective potential of the field soon after the end of inflation at <fi ~ M p is dominated by the term 
^0 4 . Oscillations of the field <p i n this theory are not sinusoidal, they are given by elliptic functions, 

but with a good accuracy one can write 4>(t) ~ $ sin(cVA where c = w 0.85. The 

Universe at that time expands as at the radiation-dominated stage: a(t) oc \fi. If one neglects the 
feedback of created 0-particles on the homogeneous field <fi(t), then its amplitude <&(£) oc a _1 (t), 
so that a<3> = const. Using a conformal time rj, exact equation for quantum fluctuations 5<j) of the 
field <p can be reduced to the Lame equation. The results remain essentially the same if we use an 
approximate equation 



d 2 (5<P k ) 



r cit At 

k 2 + 3Aa 2 $ 2 sin 2 {cVXa^r]) 5^ = , rj = / — - = — - , (4) 
an 1 1 1 J a(t) a(t) 

which leads to the Mathieu equation with A = + ^f? ~ ^Sl^ + ^-08, and q = 1.04. 

Looking at the instability chart, we see that the resonance occurs in the second band, for k 2 ~ 
3Aa 2 $ 2 . The maximal value of the coefficient //^ in this band for q ~ 1 approximately equals to 
0.07. As long as the backreaction of created particles is small, expansion of the Universe does not 
shift fluctuations away from the resonance band, and the number of produced particles grows as 
exp(2c^v / Aa$r/) ~ exp(^). 

After the time interval ~ M p " 1 A~ 1//2 | In A| , backreaction of created particles becomes significant. 
The growth of the fluctuations (0 2 ) gives rise to a contribution 3A(0 2 ) to the effective mass squared 
of the field 0, both in the equation for <p(t) and in Eq. ([|) for inhomogeneous modes. The stage 
of explosive reheating ends when (0 2 ) becomes greater than $ 2 . After that, $ 2 (0 2 ) and the 



effective frequency of oscillations is determined by the term y'3A(0 2 ). The corresponding process 

is described by Eq. (4) with A(k) = 1 + 2q + 3A jf^2\ , q = ^jj^y- In this regime q 1, and particle 
creation occurs in the narrow resonance regime in the second band with A ~ 4. Decay of the field 
in this regime is extremely slow: the amplitude $ decreases only by a factor t 1 ^ 12 faster that it 
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would decrease without any decay, due to the expansion of the Universe only, i.e., $ oc t 7 / 12 0. 
Reheating stops altogether when the presence of non-zero mass though still small as compared 

to y / 3A(0 2 ) appears enough for the expansion of the Universe to drive a mode away from the narrow 

resonance. It happens when the amplitude $ drops up to a value ~ m^/vA. 

In addition to this process, the field may decay to ^-particles. This is the leading process 
for g 2 ^> A. The equation for Xk quanta has the same form as eq. (^) with the obvious change 
A — > g 2 /3. Initially parametric resonance is broad. The values of the parameter /!& along the line 
A = 2q do not change monotonically, but typically for q ^> 1 they are 3 to 4 times greater than the 
parameter /i*. for the decay of the field into its own quanta. Therefore, this pre-heating process 
is very efficient. It ends at the moment t ~ M^A -1 / 2 ln(A/g 10 ) when $ 2 ~ (x 2 ) ~ g^^/XM 2 . The 
typical energy of created x-particles is E ~ (g 2 \) 1 ^ 4: Mp. The following evolution is essentially the 
same as that described in Sec. 3. 

5. Finally, let us consider the case with symmetry breaking. In the beginning, when the 
amplitude of oscillations is much greater than a, the theory of decay of the inflaton field is the 
same as in the case considered above. The most important part of pre-heating occurs at this 
stage. When the amplitude of the oscillations becomes smaller than m^/vA and the field begins 
oscillating near the minimum of the effective potential at <fi = a, particle production due to the 
narrow parametric resonance typically becomes very weak. The main reason for this is related to 
the backreaction of particles created at the preceding stage of pre-heating on the rate of expansion 
of the universe and on the shape of the effective potential ||. However, importance of spontaneous 
symmetry breaking for the theory of reheating should not be underestimated, since it gives rise to 
the interaction term g 2 a<px 2 which is linear in 0. Such terms are necessary for a complete decay of 
the inflaton field in accordance with the perturbation theory ([!]). 

6. In this paper we discussed the process of reheating of the universe in various inflationary 
models. We have found that decay of the inflaton field typically begins with a stage of explosive 
production of particles at a stage of a broad parametric resonance. Later the resonance becomes 
narrow, and finally this stage of decay finishes altogether. Interactions of particles produced at 
this stage, their decay into other particles and subsequent thermalization typically require much 
more time that the stage of pre-heating, since these processes are suppressed by the small values of 
coupling constants. The corresponding processes in many cases can be described by the elementary 
theory of reheating. However, this theory should be applied not to the decay of the original large 
and homogeneous oscillating inflaton field, but to the decay of particles produced at the stage of 
pre-heating, as well as to the decay of small remnants of the classical inflaton field. This makes a 
lot of difference, since typically coupling constants of interaction of the inflaton field with matter 
are extremely small, whereas coupling constants involved in the decay of other bosons can be much 
greater. As a result, the reheating temperature can be much higher than the typical temperature 
T r < 10 9 GeV which could be obtained neglecting the stage of parametric resonance ||. On the 
other hand, such processes as baryon creation after inflation occur best of all outside the state of 
thermal equilibrium. Therefore, the stage of pre-heating may play an extremely important role in 
our cosmological scenario. Another consequence of the resonance effects is an almost instantaneous 
change of equation of state from the vacuum-like one to the equation of state of relativistic matter. 
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This leads to suppression of the number of primordial black holes which could be produced after 
inflation. 
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Figure Caption: 

Fig. 1. The sketch of the stability/instability chart of the canonical Mathieu equation (3). 
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